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The staggered fermion approach to build models with chiral fermions is briefly reviewed. The method is tested 
in a U(l) model with axial vector coupling in two and four dimensions. 



1. INTRODUCTION 

An important open problem in lattice field the- 
ory is to find a satisfactory method to describe 
chiral fermions. There exist several attempts and 
proposals, for an overview see jl|, but none of 
these has so far been shown to work. Here we 
shall discuss and test an approach Q based on 
the staggered fermion method. Note that in this 
approach the notorious species doublers need not 
be decoupled, since they are used as physical spin- 
flavor degrees of freedom. 

Staggered fermion fields, denoted by Xx, do not 
carry explicit flavor and Dirac labels since these 
components are 'spread out' over the lattice. In 
the classical continuum limit one can see the usual 
Dirac and flavor structure in momentum space 
jH, and an equivalent method can also be given 
in position space. Thereto one introduces 4x4 
matrices ^"'^ defined as, 



'Xx+h, = 



b 



with 7^ = 7i ^ ■ ■ '74* and the sum running over 
the corners of a hypercube, 6^ = 0, 1. Note that 
these Vl/a; are not independent, only their Fourier 
components in the restricted momentum range 
7r/2 < < 7r/2 may be considered independent. 
The index a{K) on ^ acts like a Dirac (fiavor) in- 
dex and one can construct models involving arbi- 
trary spin-flavor couplings to Higgs or gauge fields 
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in a straightforward manner In the classical 
continuum limit (for smooth external fields and 
for the low momentum modes of the desired 
symmetry properties of the model emerge. The 
important issue is if the same can be achieved 
when the effect of quantum fluctuations is taken 
into account. 

This approach has been successfully applied to 
a study of a fermion-Higgs model Here we 
consider a U(l) model with axial- vector coupling 
in two and four dimensions, defined by the action 
in the ^E* notation, {d — 4), 

= -T. [^x1^{U^xPl + U;,Pn)^,+^ 

- '^x+i.7^(U;,PL + U^xPB)'fx]. (2) 

By working out the trace in (j^) one finds the 
action for the staggered field. 



32 / y 



/ , '^x-bXxXx+fi 
b 



h.C.) 



{iVbxVlix+n E ^x-bXxXx+c-b+i. + h.C.) 
b-\-c—n 



(3) 



where ctj: = Ret/^j, and = ImU^ix- The sign 
factor ry^a: = (— 1)^^"^ \-x^-i represents the Dirac 

7a., V5x = »7lx»72x+i---^4x+i + ...+4 ^ (-1)^1+^^ 

represents 75 and n — (1,1,1,1). In the classical 
continuum limit this action describes four (two) 
flavors of axially coupled Dirac fermions in four 
(two) dimensions. 
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In the 4" form the action appears to be gauge 
invariant under [/ ^ J7^, ^' with 

f/^. = 5.C/p.5^+A, (4) 

n = {g.PL + g^PR)^.. (5) 

However, since the ^'a; are not independent, such 
a local transformation cannot be carried over to 
the X field, and the action is not gauge invariant. 

2. GAUGE (NON)INVARIANCE 

Even though the action (||) is not gauge in- 
variant, it is well-known ||^ that the model is 
closely related to a gauge invariant model. Let 
S = Sf + Set be the gauge non-invariant action, 
consisting of the term (^) plus additional coun- 
terterms, denoted by Set- The partition function 
can be written as, 

^-i""' I Dx'Dx'DVDUe'^''''''^'''^'\ (6) 

Here we have used invariance of the Haar measure 
DU = DU\ g = V\ and J DV = 1. By explic- 
itly writing the 'longitudinal gauge mode' V, the 
action S{U^ , Xj x) is seen to be gauge invariant 
under (||) and 14 —^gxVx- 

In the scaling region this gauge mode V may 
decouple from the physical states (if the fermion 
content is anomaly free), which is to be achieved 
by choosing the appropriate counterterms in Set- 
In particular we add a mass term for the gauge 
field. 

Set = Kj2K^ + h.c.), (7) 
= t,J2iV:U^=cV,+f, + h.c.). (8) 

The V field can also be viewed as a radially frozen 
Higgs field. Then acts like a kinetic term 
for the Higgs field. One expects a critical line 
K — Kc{P) (with l3 the bare gauge coupling) and 
for K Kc the Higgs mode will be present as a 
physical state. By keeping k away from Kc in the 
symmetric phase, the scalar mass will be of the 
order of the cut-off and V will decouple. 



The attractive feature here is that gauge fix- 
ing is avoided. Alternatively one might try to fix 
the gauge and then remove undesired couplings 
of the V field by adding suitable counterterms, 
as in the 'Rome approach' One then has to 
worry about technical and Gribov problems with 
non-perturbative gauge fixing. 

Using y as a Higgs field, we may think of U' as 
the gauge field in the unitary gauge, and similar 
for 4"' and x'- It is natural to add a fermion mass 
term 

5y = -y^Tr*l*^, (9) 

X 

which would turn into a Yukawa couphng 
-vEx^^^AV^'Pr + {V*)'Pl)^. if % = 
{V*Pl + VxPR)'i'x could be implemented as a 
transformation of variables. 

3. EXTERNAL GAUGE FIELDS 

We shall first test our model for smooth ex- 
ternal gauge fields, i.e. U ~ 1 everywhere, and 

= 1, in two dimensions. Since our model has 
axial-vector coupling, the vector current has an 
anomaly. One therefore expects the relation, 

id'^J^^o.)^ = ^Fx, (10) 

where the notation (•)^ means integration over 
the X field. The vector current (given in ^E" 
notation for simplicity) and field strength F are 
defined as, 

J^, = iTii^xlf.U^x'ifx+f. + h.c.), (11) 
Fx+n/2 = -ilogU^xU^^^iUl^^U;,, (12) 

As an example we consider a spatially constant 
U field defined by, 

Uix = exp{iAi sin(27rt/i)), U2x = 1. (13) 

Lattice sites are denoted by a; = {z,t) and the 
lattice size is L^. Fig. 1 shows the divergence of 
the current (the left hand side of eq. (|lO|), denoted 
by □) and the anomaly (full line). Clearly the 
correspondence is excellent and this remains the 
case for amplitudes Ai^O.2. For larger values one 
begins to see deviations. 



face effect. 
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Figure 1. Current divergence (□) and anomaly 
(full line) for the smooth gauge field ([T^), L — 24. 
The (*) includes averaging over V field fiuctua- 
tions with k = 0.6, L = 8. 



Our second test case uses an external gauge 
field with non-trivial topology. For the field (13) 
the topological charge Q = Fx/2t: is zero. We 
can construct a gauge field with non-zero Q and 
constant field strength Fx = by writing 



[/la: = exp(iFt), U2x ^ I ,t ^ 1,- ■ ■ ,L - I; 

U2x — exp{iFLz), t ~ L. 



(14) 



For small Q/L, the link field U is close to one ev- 
erywhere except on the time slice t = L which con- 
tains a transition function. This transition func- 
tion makes U 'rough' (discontinuous) &i t = L. 
With a gauge invariant model such a transition 
function is invisible, but since our action and the 
current ( p]) lack gauge invariance we may expect 
deviations from the anomaly equation (10) at this 
timeslice. This is illustrated in fig. 2, which is the 
same as fig. 1 but now for the gauge field (p^. 
The presence of the slice with the transition func- 
tion is clearly visible and its disturbance extends 
over roughly 6 time slices. We have checked that 
on a larger lattice the disturbance remains con- 
fined to this number of time slices and one might 
argue that in the infinite volume limit the effect of 
the transition function becomes an invisible sur- 



Figure 2. Same as fig. 1, now for the gauge field 
([l4| ) with Q = 1 and without V field fluctuations. 



These results show that for smooth fields the 
correct anomaly structure is reproduced by the 
staggered fermions. Also non-trivial topological 
sectors can presumably be incorporated in the 
proper way, but the strong effect of the rough- 
ness produced by the transition functions hints 
at difficulties for fields subject to quantum fluc- 
tuations. 

4. FLUCTUATING V FIELD 

Suppose we decompose the gauge field in a 
'transverse' and a 'longitudinal' part, U'^^ = 
VxUf^x^*^^, where C/*'' satisfies the (lattice) Lan- 
dau gauge condition. The fluctuations of [/*'' 
are governed by a gauge coupling /3 and get sup- 
pressed for large /3. If not used as a physical Higgs 
field, the V modes should be unobservable in the 
low energy physics, if the model becomes gauge 
invariant. This can be tested already for [/*'' 1. 

We have investigated the resulting fermion- 
scalar model for y ^ 0. For y ^ and k 
strictly smaller than Kc the V field should de- 
couple. In the quenched approximation wc have 
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computed the fermion propagators {^x'^y) and 

where = (V^Pr + V^Pl)^. Only 
the first propagator should then show a physi- 
cal particle pole (in the symmetric phase). We 
found encouraging results in the broken phase but 
ambiguous results, hampered by huge statistical 
fluctuations, in the symmetric phase. 

As an alternative test we add the longitudinal 
modes to the smooth external gauge field jl^), 
Uf^x U'^x = VxUf_,xV*^^, and compute the 
anomalous divergence relation (|l^) with the path 
integral average over V included {y — 0). For 
large k^I we then reproduce the result for the 
fixed background, up to statistical errors. For 
smaller k, approaching the phase transition at 
Kc ~ 0.5 we find that the current gets renormal- 
ized. The result in fig. 1 (denoted by *) shows 
that the current strength is reduced by a factor 
« 0.7 at K = 0.6. Such a renormalization is to 
be expected for staggered currents and has been 
studied previously in QCD, see e.g. ||]. 

For smaller values of n, in particular inside the 
symmetric phase, the magnitude of the statisti- 
cal fluctuations grows dramatically and prohibits 
computation of the average current divergence. 

5. DYNAMICAL FERMIONS 

Since the back reaction of the fermions on the V 
field should reduce the troublesome fiuctuations, 
we have also studied the scalar-fermion model in 
four dimensions with dynamical fermions. Here 
we focus on the phase diagram of the model in 
the y — K plane. For large y one can use a hopping 
expansion to see that the fermion field ^' in this 
unitary gauge has a pole whereas the fermion field 
* does not. If this persists also for small y the 
model fails. If the model is viable we must see a 
cross over at some (small) value of y. 

Unfortunately the problems which hampered 
the quenched computations appear to haunt us 
here as well: For small values of y'^OA the model 
has extremely long living metastable states, rem- 
iniscent of spin glass dynamics and we have so 
far not been able to compute the phase structure 
in the region y^O.4. The phase diagram must 
be symmetric for k — > — k, which we can use to 
check our results. In fig. 3 we show the measured 



Figure 3. Possible phase diagrams in the k — 
y plane of the four dimensional axial U(l) 
model. PM (FM) indicates the symmetric (bro- 
ken) phase, the PM phase at small y is absent if 
the solid lines are true. 



phase transitions for large y and two speculations 
for the structure at small y. One possibility, the 
solid lines, represent first order transitions and 
would rule out our model. A second possibility is 
indicated by the dashed, second order lines and 
in this scenario our model is still alive. 

The numerical calculations were performed on the 
CRAY Y-MP4/464 at SARA, Amsterdam. This 
research was supported by the "Stichting voor 
Fundamenteel Onderzoek der Materie (FOM)" 
and by the "Stichting Nationale Computer Fa- 
ciliteiten (NCF)" . 
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